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ABSTRACT 


The Hartree-Fock method is used to study the even- even N Z nuclei in the 2s- Id 
shell. Basis functions with correct asymptotic behavior are generated by use of a Wood- 
Saxon well. The results are compared with calculations using the conventional harmonic 
oscillator basis. It is found that most of the Wood-Saxon results may be obtained with a 
harmonic oscillator basis if the oscillator parameter is carefully chosen. 



HARTREE-FOCK CALCULATIONS WITH WOOD-SAXON BASIS FUNCTIONS 
by William F. Ford and Richard C. Braley 
Lewis Research Center 

SUMMARY 

In order to determine the sensitivity of various nuclear properties to the asymptotic 
part of the nuclear wave function, calculations have been made using both the conven- 
tional harmonic oscillator single-particle functions and the more realistic Wood- Saxon 
functions. Hartree-Fock theory was used to obtain nuclear wave functions for the five 
even- even N = Z nuclei in the 2s- Id shell. Intrinsic nuclear energies, radii, and quad- 
ruple moments were calculated, as well as energy gaps and single-particle densities. 
Considerable differences between harmonic oscillator and Wood-Saxon results were 
found when the usual value (v = 0. 35 fm ) for the oscillator parameter was used; how- 
ever, close agreement could be obtained in most of the properties studied by using a 
different value (u = 0. 27 fm"^). 


INTRODUCTION 

In recent years the Hartree-Fock (HF) method has received considerable attention 
in studies of the energy levels of deformed nuclei (ref. 1). This method is essentially a 
variational procedure, based on minimization of the energy of a many-particle system. 
While such a procedure may provide a sensitive test of the wave function in the nuclear 
interior, it is expected that the shape of the wave function beyond the nuclear surface 
will not affect nuclear energies significantly. For this reason, HF calculations up until 
now have made use of single-particle harmonic oscillator basis functions, in spite of 
their obviously incorrect asymptotic behavior. 

There are, however, many quantities of interest which are expected to be much more 
sensitive than the total energy to the asymptotic part of the wave function. Electro- 
magnetic transition rates, multipole moments, and form factors for elastic and inelastic 
scattering are examples which come readily to mind. It does not seem likely that very 
accurate calculations of such quantities can be performed using HF wave functions with a 
harmonic oscillator basis. Thus, a HF calculation using the basis functions of che more 


realistic Wood-Saxon potential is of interest both for the energy spectrum it yields and 
for the effect it shows of the long-range part of the wave function on other nuclear prop- 
erties. 

In this report are presented the results of a considerably restricted version of such 
a program. The HF variational calculation is carried out for five nuclei in the 2s- Id 
shell by using both Wood-Saxon (WS) and harmonic oscillator (HO) basis functions. A 
comparison is made of such properties as orbital energies, energy gap, and total (HF) 
energy. In addition, nuclear densities, radii, and intrinsic quadrupole moments are 
calculated and compared. 


SYMBOLS 



e x 

h 


J 

(jm) 


(^mp j 2 m 2 | JM> 


or 


li> 


fi. 


s a Ja 


f b b J b 
L S J 


(m) 


P ij 

r ij 


diffuseness of Wood-Saxon well 

expansion coefficients for the Hartree-Fock orbits 

Hartree-Fock energy 

Hartree-Fock single-particle energy 
Hartree-Fock Hamiltonian 
independent-particle Hamiltonian 
V2J + 1 

single -particle total angular momentum and its projection onto 
z-axis 

Clebsch-Gordan coefficient 

reduced matrix element 

single -particle orbital angular momentum 


symmetric 9-j coefficient 
n! 

m!(n - m)I 

two-particle exchange operator 
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total spin angular momentum and projection onto z-axis 


(SMg) 

U(jij 2 j 4 33’ J 12 J ) U-coeffieient of Jahn (see ref. 9); defined in terms of Racah coef- 

ficient as [(2J 12 + 1)(2J + l)] 1 / 2 W(J 1 J 2 i 4 j 3 ;J 12 J) 

“a<‘> ’■a< r l> T r„ a V2J„ (1)X l/2 (l) 

a a 

u„(r.) radial part of basis function 

a l 

V(ij) residual two-body force 

V Q strength of residual two-body force and/or Wood-Saxon well 

v(i) nonlocal part of the Hartree-Fock Hamiltonian 

Racah coefficient 

Y^\i) spherical harmonic 


Y ?s jW 

A 

M 

v 

p( r) 

PyM 

a(i) 

T a 

r dr. 
^ J 

$ 

X SM g 

Xi/ 2 (i) 


E m-j m 

(m 1 ,sm 2 |jm) (i)X g 


2 (i) 


m l m 2 

energy gap between occupied and unoccupied orbits 
nucleon-nucleon force range 
Harmonic oscillator parameter 

single-particle density function which includes core effects 
single-particle density function for valence particles only 
Pauli spin operator for particle i 
z-projection of single-particle isotopic spin 


in 

integral over all coordinates of j particle 

determinantal wave function for A-particle system 
Hartree-Fock orbital function of i cn particle 
total spin wave function for two-particle system 

1L 

spin wave function for i Ln particle 

1 1 

isospin function for i particle 
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THEORY 


In the HF method, the nuclear wave function is approximated by the Slater deter - 
minant 


cE> = (A!)" 1/2 det{ n (i)} (1) 

whose orbital functions cp x ( i) are determined by minimizing the expectation value of the 
nuclear Hamiltonian. This variational problem leads to the eigenvalue equation 

hcp x = (h Q + v)<p x = e x <p x (2) 

where 

v(i) =- f dT^*(j)V(ij)(l - P tj )«p M (j) (3) 

2 H 

The operator h Q is usually taken to be the Hamiltonian for a single nucleon moving in the 
average nuclear field, and V represents the residual two-body force. The presence 
of the exchange operator is a natural consequence of the antisymmetry of the deter- 
minant. The summation in equation (3) is taken over all the orbits which are to be varied. 

Since equations (2) and (3) are too difficult for a direct solution, approximation 
methods are required. Usually one begins by selecting a suitable set of basis functions 
and then making the expansion 


■?*<» = £ C aV» («> 

a 

The problem then reduces to diagonalizing h in the space spanned by the basis functions 

u . In practice, this is accomplished in an approximate way by first truncating the space, 

then guessing at an initial set C , and finally iterating until two successive diagonaliza- 

A ^ 

tions yield the same set C„. (Often "a” will be written for the set j m , the set j , 

a da a. 

or even the set n i j when no confusion will result. ) 

3. d d 

In order to carry out this procedure, one needs matrix elements of h. It is gener- 
ally assumed that h Q is diagonal with eigenvalues deduced from experimental spectra. 
The matrix elements of v may be expressed (ref. 1) in terms of the coupled two-body 
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matrix elements of V: 


<a. | v | c> =i ZE(cfc;E(] A 

’ j b m bl JM)(j c m c’ j d m dl JM> 

* li bd JT 


V 5 Tb|TMr )(l T dl™r)< ab l v l' rd > J T 


( 5 ) 


(The tilde indicates that the product function u (i)u,(j) has been antisymmetrized.) The 

cd) JT are of course independent of the orbital functions, 
and they are therefore ideally suited for use as input in a general program designed to 
solve equation (2). 

Calculation of the two-body matrix elements presents some problems which stem 
from the fact that V(ij) is a function of the separation between nucleons i and j. If 
HO functions are used, one can overcome the difficulty by transforming to center-of- 
mass and relative coordinates, and the resulting integrals can be expressed in terms of 
Moshinsky brackets (ref. 2). In the general case, however, such a transformation is 
not possible, and the most promising approach seems to be an expansion of V(ij). 

For the moment the isotopic spin dependence is ignored. The coupled matrix ele- 
ment may then be written 


<ab| V|cd>j = 1 2 > I v 0 2 >l 1 2> " ^M 121 ^ 


( 6 ) 


where 


f>JM< 12 > = X <ia m a.!b m bl JM > u a< 1 > u b< 2 > = (- l) la+J 2 1) 

m a m b 

After switching from j - j to L - S coupling (ref. 3) by means of 

«>JM = J] <LM L’ SM sI JM> ' I ’LM l %M s 

m l m s 


( 7 ) 


we obtain 
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(8) 


(^ 12 >i v < 12 >k$S!i< 1 ») = Z <sj^j 

X (^LM l ( 12 )x SM s (12 ’ I v<12) l' X Lrvl L (lJ,X SM s flJ, \ 

where M-^ + Mg = M and 

pa >/ 2 C 

<sj s £i y b z b 1x2 J < 9 > 

L S J 

J 

In this derivation V is assumed invariant under space and spin rotations separately. 
Usually a spin dependence of the form 

1 

V(12) = ^ V s ,(12)5g,(l) • ct s ,(2) (10) 

S’=0 

is also assumed where ov. is the unit operator in spin space and <71 is the ordinary 

u 1 4-9 

Pauli spin operator. Then, since Xg(12) = (-1) Xg(21), 

(ab| V|cd) j = QlsjQlSJ X/ A SS’ 

LS S’ 

X I Vg,(12) | + (-1) 

where 

A ss? = <S| |5g,( 1) • 5g,(2)||s) 

= 2(2S- + D(-l> S ' l w(iiii|s's) = 1 - 45 S0 5 S ,J (12) 

nh 

Next the spherical harmonics in <p LM and (Pt\^ are decoupled and then recoupled 
according to 
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(13a) 




a a c c 


L'M’ 


Y ‘h m h (2)Y *H m H (2) = ^ <Z dl l V L” 

b b d d L m M m 


l i c><'c m c' L ' M, Ua“a> Y L'M’< 1 ) 

| i b > <i b m b , L"M" | ; d m d ) Y l „ m „( 2 ) 


(13b) 


Now, after a little more Racah algebra, the coupled matrix element may be written 

(abj V) cd) j = ^ Q LSJ Q LSJ X A SS’ ^LL^I/S^ 

LS L’S’ 

+ (-l) L+S+Zc+ ^ d B^, cb F^ cb ] (14) 

where 

B & db s « a l l Y L'l l'c>< ! bl l y L’l I ' b > u < J c I/Li bl Vd> (15) 

F S’ db = (l»LM( 1 )l V a< 12 )l« , LM< 2 >) < 16 > 

and 


~ u a^ r l^ u c^ r P Y LM^ r l^ ( 17 ) 

In the derivation of equation (14) Vg is assumed central so that the sums over L" and 
M" reduce to a single term. 

The final step is the inclusion of isotopic spin dependence. It is assumed that 

1 

V gt (12) = X V s , T ,(12)r T ,(l) • t t ,(2) (18) 

T'=0 

where Tq and Tj are defined like <?q and but in isotopic spin space. The treat- 
ment is evidently exactly like that given for the spin dependence. With the customary 
assumption that 
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Vs t T i( 1 2 ) - ^S'T ,V ^ l r l - r 2 |) 


(19) 


where /3g, T , measures the relative strengths of V gf T , , the final result is obtained: 


(ab j V | cd) ^ - ^Qlsj Q LSJ 23 A SS' A TT’%T' 
LS L'S’T’ 


R ac, db„ac, db 
B LL* F L’ 


- (- 1 ) 


L+S+T+Z c +Z 



( 20 ) 


where 


fL C ’ db = (21) 

The method used to evaluate integral (21) is based on a technique due to Sawaguri 
and Tobocman (ref. 4), who developed a general formalism for evaluation of six- 
dimensional integrals. A simplified version adequate for the work herein is obtained as 
follows . Let 


/imn® = 

be a complete set of functions in terms of which the expansion 

V(|r- r 2 |) = \n 2 lhmn^n^2) 
l m 
n l n 2 


is made. 


The desired integral (21) may then be written 

oo 


F 


ac, db 
L 


z 

n,n„ 


V T? aC 

V n 1 n 2 L F n 1 L 



where 


( 22 ) 


(23) 


(24) 
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(25) 


Kl ’ X" u*( r ^(r)u c (r)r 2 dr 


The usefulness of the method depends on whether the series in equation (24) is 


rapidly convergent and whether values of V. 
Tobocman showed that for the choice 


n iV 


are easy to obtain. Sawaguri and 


and 


2„2_2, 


^(r) =& l n (a,pr) = _L (/3r) Z e"“ ^ r /^F^-n ;l + 3/2;/3 2 r 2 ) 


(26) 


one gets 


n l n 2 


V = \ ' \ V 

n lV m l n l m 2 n 2 m l 


+m2+Z 


m ^=0 m 2=0 


(27) 


where 


V = 
m 


\V2 - a/ (-of) m «/) ' 


i^Hr 2 dr 


(28) 


mn 


i^r /2 o 


r n +1 + - 


r/z + 2 \r/m + 1 +- 


(29) 


In practice, for appropriate choices of a and /3, the series in equation (24) can be made 
to converge sufficiently fast to make the entire scheme an excellent computational device. 


RESULTS 

The theory just described has been applied to the five s-d shell nuclei neon 20, mag- 
nesium 24, silicon 28, sulfur 32, and argon 36 (Ne 2 ®, Mg 2 ^, Si 2 ®, S 22 , and A 2 ®). The 
HF wave functions are assumed to have axial symmetry (ref. 5), and they are formed by 
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putting particles into deformed orbits about a spherical, closed-shell, inert oxygen 16 

1 d 

(O ) core. Each deformed orbit is described as a mixture of ld^* 2 s i/ 2> and ld 3/2 
basis vectors, each with the appropriate spin projection m^. The coefficients for the 
negative m^ orbits are related to those for the positive m^ orbits by time reversal 
symmetry: 


C*(-m x ) = cjV x )(-l) 3 ' mX (30) 

The eigenvalues of hr used in the calculation are taken from the experimental spec- 
17 ° 

trum of oxygen 17 (O ), and they are listed in the fourth column of table I. The residual 
two-body force V(ij) is a Rosenf eld mixture with Gaussian shape, much like that used by 
Ripka (ref. 1): 


V(ij) = V Q [o. 3 + 0.7 ffjd) • aj(j)] e ^ (31) 

3 

with V - -50 MeV (instead of Ripka's V = -55 MeV) and ju = 0. 29 per square femtom- 
eter (fm - ^). 

For a discussion of the basis vectors themselves, two sets of harmonic oscillator 
basis functions were used and one set of Wood-Saxon basis functions (ref. 6). The latter 
were chosen by varying the parameters V Q and a in the Wood-Saxon potential 



so as to reproduce as closely as possible the experimental eigenvalues of h Q as given 

in table I. The eigenvalues thus obtained are given in the third column of table I. 

Choosing the harmonic oscillator basis functions poses something of a problem - 

namely, what criteria to use. At length it was decided that one set of functions (HO - 1) 

should be chosen on the basis of similarity to the Wood-Saxon basis functions within the 

nuclear interior, as measured by comparison of eigenvalues and location of peaks and 

_2 

nodes in the functions. For this set an optimum choice was v = 0. 35 fm . This choice 
v is based on an approximate method for calculating the "correct" nuclear radius 
(ref. 7). 
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TABLE I. - OXYGEN 17 SINGLE- PARTICLE ENERGIES 


Orbit 

Energy, MeV 


Harmonic oscillator parameter, 

, -2 
v, fm 

Wood- Saxon 

Experiment 


v = 0. 27 fm* 2 

v = 0. 35 fm' 2 



ls l/2 

-26. 2 

-32. 6 

-34. 68 

— 35 

^3/2 

-15.0 

-18.2 

-20. 21 

-21. 83 

^1/2 

-15.0 

-18. 2 

-16. 53 

-15. 67 

ld 5/2 

-3.7 

-3. 7 

-6.02 

-4. 14 

2s 1/2 

-3.7 

-3. 7 

-4. 54 

-3. 27 

ld 3/2 

-3.7 

-3. 7 

i 

o 

cn 

CO 

0. 94 


A second set (HO - 2) was chosen eventually on the basis of similarity to the Wood- 
Saxon basis functions as measured by comparison of the resulting nuclear binding 
energies, orbital energies, energy gaps, and intrinsic quadrupole moments. For this 
set, v = 0. 27 fm -2 . 

In table I the energy eigenvalues of the various basis functions are compared, while 
in figures 1 to 5 their shapes are compared. Little can be determined from the energy 
eigenvalues, but figures 1 to 5 reveal clearly for the ls-ip shell functions that HO-1 is 
closer to WS than HO-2 is. In the 2s- Id shell the situation is somewhat ambiguous: 

HO - 1 has the same peaks and nodes as WS, as before, whereas HO-2 has different peaks 
and nodes but has perhaps an overall better match in magnitude. 

The results of HF calculations using the three sets of basis functions just described 
are compared in tables n(a) and (b). On the basis of the binding energies and the gaps, 
it is evident that HO-2 gives a better match to the WS results than HO-1. 


Root-Mean-Square Radii 

The rms radii are obtained from evaluation of the expectation value of the operator 


R 2 



i=l 


A 


(34) 
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Parameter 


Wood Saxon 

Harmonic oscillator, v = 0. 35 

Harmonic oscillator, v * 0. 27 



~N 

\ 



3 4 5 6 

Radius, r, fm (or F) 


(el 

Figure 1. - Concluded. 








TABLE II. - ENERGIES AND WAVE FUNC TIONS FOR AXIALLY SYMMETRIC, EVEN- EVEN, 2s- Id SHELL NUCLEI 
[Results using a Wood-Saxon basis are in boldface type; light type corresponds to a harmonic oscillator basis with a stated harmonic oscillator parameter 14] 


(a) Harmonic oscillator parameter, v - 0. 35 fm' 


Nucleus 

Hartree- Fock 
energy. 

e hf 

MeV 

Energy gap. 
A, 

MeV 

Single- particle 
Hartree- Fock 
energy, 

e x> 

MeV 

Projection of 
angular momentum 
on z-axis, 

m x 

Projection of the Hartree- Fock orbits on spherical basis, 

c ? 

ld 5/2, 1/2 

ld 5/2, 3/2 

ld 5/2, 1/2 

2s l/2, 1/2 

ld 3/2, 1/2 

ld 3/2, 3/2 

20 

Ne 

-40.21 

10.6 

-17.01 

1/2 

0.6781 



-0.6342 

-0.3715 



-35.27 

8.5 

-14.46 

1/2 

.7782 



-. 4943 

-.3874 


Mg 24 

-72.28 

4.9 

-15.88 

3/2 


0.9559 




-0.2938 


-68.26 

5.0 

-14.78 

3/2 


.9550 




-. 2835 




-13.83 

1/2 

0.8414 



-0.1080 

0.5294 





-12. 78 

1/2 

. 8552 



-.0997 

. 5086 


Si 28 

-132.45 

8.2 

-21.93 

1/2 

0.4798 



0.8320 

-0.2787 



-116.80 

7. 1 

-18. 37 

5/2 



1.000 







-19.95 

5/2 



1.000 







-17. 16 

1/2 

0. 6066 



0. 7345 

-0. 3040 





-15.86 

3/2 


0.6440 




0.7650 




-13.98 

3/2 


. 6959 




. 7182 

s 32 

-179.64 

5.9 

-22.66 

1/2 

0.4083 



0.8830 

-0.2315 



-157.89 

5.0 

-18.35 

5/2 



1.000 







-20.24 

5/2 



1.000 







-18. 11 

1/2 

0. 8167 



0. 4404 

-0. 3730 





-18.55 

1/2 

0.7965 



-0.4685 ■ 

-0.3823 





-15.14 

1/2 

-.4162 



.8972 

. 1481 





-16.22 

3/2 


0.7030 




0.7112 




-14.35 

3/2 


.7889 




.6145 

Ar 36 

-238.56 

8.1 

-25.07 

1/2 

0.3607 



0.9268 

-0.1040 



-209. 66 

6.9 

-21.41 

5/2 



1.000 







-23.05 

5/2 



1.000 







-20.38 

3/2 


0.9968 




0. 0801 




-22.27 

3/2 


0.9969 




0.0781 




-19.47 

1/2 

0. 6691 



0. 7419 

0.0431 





-18.94 

1/2 

0.7411 



-0.2172 

0.6352 





-16.16 

1/2 

. 5498 



-. 5332 

. 6429 





-17.18 

3/2 


-0.0781 




0.9969 




-14. 57 

3/2 


-.0801 




.9968 


(b) Harmonic oscillator parameter, v = 0. 27 fm 


-2 


Ne 20 

- 34.82 

-35. 27 

8.2 

8. 5 

- 14.27 

-14.46 

1/2 

1/2 

0.7273 

.7782 



- 0.5729 

-.4943 

- 0.3780 

-. 3874 


Ms 24 

- 65.99 

4.8 

- 14.14 

3/2 


0.9562 




- 0.2928 

IVlg 

-68. 26 

5.0 

-14.78 

3/2 


.9550 




-. 2835 




- 12.52 

1/2 

0.8333 



- 0.0809 

0.5469 





-12.78 

1/2 

. 8552 



-.0997 

. 5086 


Si 28 

- 116.41 

7.0 

- 17.94 

1/2 

0.5290 



0.7967 

- 0.2922 



-116. 80 

7. 1 

-18. 37 

5/2 



1.000 







- 17.55 

5/2 



1.000 







-17. 16 

1/2 

0. 6066 



0.7345 

-0. 3040 





- 14.07 

3/2 


0.6682 




0.7440 




-13.98 

3/2 


. 6959 




. 7182 

s 32 

- 159.06 

5.0 

- 18.76 

1/2 

0.5757 



0.7598 

- 0.3020 



-157. 89 

5.0 

-18. 35 

5/2 



1.000 







- 17.76 

5/2 



1.000 







-18. 11 

1/2 

0. 8167 



0. 4404 

-0. 3730 





- 16.11 

1/2 

0.6928 



- 0.6495 

- 0.3133 





-15. 14 

1/2 

-. 4162 



. 8972 

. 1481 





- 14.48 

3/2 


0.7333 




0.6799 




-14. 35 

3/2 


. 7889 




. 6145 

Ar 36 

- 211.89 

7.1 

- 20.68 

5/2 



1.000 





-209. 66 

6. 9 

-21.41 

5/2 



1.000 







- 20.40 

1/2 

0.4868 



0.8716 

- 0.0576 





-20.38 

3/2 


0.9968 




0. 0801 




- 19.81 

3/2 


0.9960 




0.0895 




-19.47 

1/2 

0. 6691 



0. 7419 

0.0431 





- 16.87 

1/2 

0.6828 



- 0.3386 

0.6475 





-16. 16 

1/2 

. 5498 



-. 5332 

.6429 





- 15.42 

3/2 


- 0.0895 




0.9960 




-14. 57 

3/2 


-.0801 




. 9968 


between HF states <3>. The results, shown in figure 2, differ considerably for the three 
sets of basis functions, and only the WS results bear a reasonable resemblance to exper- 
iment and then only in an average sense. As the mass number increases toward 40, the 
WS results begin to deviate even more from experiment. The HO results are quite poor 
in both cases. However, the comparison of all these results with experiment would be 
much more meaningful if the expectation values were taken with respect to states of good 
angular momentum. The results are, nonetheless, significant for nuclear matter studies 


Intrinsic Quadrupole Moment 


The intrinsic quadrupole moment is defined as the expectation value of the operator 


Qr 



(35) 


with respect to HF states $. The operator is an isospin operator with eigenvalues 
+1 and -1 for protons and neutrons, respectively; e is the proton charge. In table HI 
the results are given for the intrinsic quadrupole moments. The set HO-2 yields a much 
better match to the WS than does HO-1. The sign of (Q 0 ) for any particular nucleus 
does not depend on the choice of basis, and the signs are in agreement with the accepted 
type of deformation (i. e. , prolate or oblate) for these nuclei. 


TABLE HI. - INTRINSIC QUADRUPOLE MOMENT 


Nucleus 

Quadrupole moment, 
<Q o >, e fm 2 



Harmonic oscillator parameter, 

* -2 
v, fm 

Wood- Saxon 


v - 0. 27 fm ^ 

v = 0. 35 fm -2 


Ne 20 

29.00 

22. 30 

29. 24 

Mg 24 

15.80 

12. 60 

14.90 

Si 28 

-42. 59 

-33. 18 

-42. 51 

g 32 

-12.74 

-10.25 

-11.74 

A 36 

-25.92 

-19. 17 

-29. 65 
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Nuclear Density 


The nuclear density function is defined as the expectation value of 

A 

P = 5(r - r 4 ) (36) 

i=l 

Since the HF calculations of the type presented here test only the extra- core wave func- 
tions, two types of density functions are shown for each of the nuclei. The valence den- 
sity p y is the expectation value of p with respect to wave functions which do not spec- 
ifically include the spherical core functions. The total density p(r) does contain the 
core effects. In figures 3 to 7 the density functions are presented. In each of the figures, 
parts (a) and (b) refer to valence and total density, respectively. The arrows which 
appear on each of the figures indicate the position of 1.25 A^'^, where A is the mass 
number. Examination of figures 3 to 7 reveals that both sets of density functions (with 
the exception of Mg^) show consistent behavior as one goes from Ne^® to A^®. For the 
moment the problem of Mg is set aside, and the general trends for the valence and 
total densities for the other nuclei are considered. 



(a) Valence density. (b) Total density. 

Figure 6. - Densities for sulfur 32. 
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Valence densities. - It is seen from figures 3(a), 5(a), 6(a), and 7(a) that for 
o 

0 ^ r ^ 1. 5 the HO-2 (v - 0. 27 fm ) results show better agreement with the WS case 
than does the HO-1 (v - 0. 35 fm - ^) result. However, for 1. 5 ^ r ^ 4. 2, HO-1 agrees 
more closely with WS. In the region beyond about 4. 2 fm, the HO-1 densities fall off 
more rapidly than both the WS and HO-2 densities. The three curves focus at about 
3. 3 fm for all four nuclei. 

Total densities . - The total densities are presented in figures 3(b), 5(b), 6(b), and 

7(b). For the WS case the flattening near the nuclear center is washed out as A in 

increases, and all of these total density functions have the same general shape. In the 

o 

region 1. 5 ^ r ^ 4. 3, the HO-1 (v = 0. 35 fm ) results match the WS results quite well, 
but drop off rapidly beyond r « 4. 5. Beyond r = 5 fm it is observed that the HO-2 
( v = 0. 27 fm ) results agree more closely with WS as might be expected since this is the 
same region in which the HO-2 wave functions agree most closely with those of the WS 
potential. 

24 

For the case of Mg , an examination of figure 4(a) perhaps gives an indication as to 
why Mg^ should be treated separately. Certainly the valence density has a peculiar ap- 
pearance in view of the consistent behavior of the other densities. However, Bar-Touv 

24 

and Kelson (ref. 8) performed calculations which indicate that Mg is not axially sym- 



Figure 7. - Densities for argon 36. 
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metric, so Mg might be expected to have properties which differ from those of axially 

QO 

symmetric nuclei. Based on the findings of Bar-Touv and Kelson, S would also be ex- 

OO 

pected to behave similarly, since the same set of calculations indicated that S is like- 
wise asymmetric. The size of the gaps in the HF spectra (see tables H(a)and (b)) indicate 

32 24 

that the assumption of axial symmetry is equally poor for S and Mg . That is, since 
the energy gap between occupied and unoccupied orbits is quite a bit smaller for these 
nuclei than for the others, a perturbation theory calculation would yield larger effects of 
axial asymmetry for these nuclei than for the others. 

24 

In addition, one of the results of attempting to treat Mg as axially symmetric is 
that the coefficient of the 2s^^2 basis vector turns out to be abnormally small, ~0. 10. 

Consequently, the valence density is drastically reduced at the origin. This does not 

32 24 

occur for S , nor does it occur in an asymmetric treatment (ref. 8) of Mg . 

If asymmetry is assumed, calculations of the energy spectra and other properties of 

these nuclei, should be more informative than the results just given. However, the 

24 

previous results do substantiate the belief that Mg has to be treated differently than the 
other 2s- Id shell nuclei. 


SUMMARY OF RESULTS 

The intrinsic energy spectra, radii, and quadrupole moments have been studied for 
the five even-even nuclei in the 2s- Id shell. The nuclear wave functions were obtained 
using the HF method. The HF calculations were made using the Wood- Saxon basis func- 
tions as well as the more conventional harmonic oscillator basis. It was found that the 
asymptotic behavior or the WS basis functions could be simulated with harmonic oscillator 
functions if the oscillator constant were properly chosen. Such oscillator functions were 
also found to yield the best match to the Woods-Saxon HF spectra, indicating an important 

dependence on the asymptotic part of the wave function. As expected, Mg (and to a 
32 

lesser extent S ) was found to have rather peculiar properties when assumed to be axially 
symmetric. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 25, 1968, 

129-02-07-07-22. 
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